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Abstract 



A quantum particle interacting with a thin solenoid and a magnetic flux 
is described by a five-parameter family of Hamilton operators, obtained 
via the method of self-adjoint extensions. One of the parameters, the 
value of the flux, corresponds to the Aharonov-Bohm effect; the other 
four parameters correspond to the strength of a singular potential bar- 
rier. The spectrum and eigenstates are computed and the scattering 
problem is solved. 
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Introduction 



The purpose of this paper is to obtain and study the most general family of oper- 
ators which describe the essential features of a quantum mechanical particle under 
the joint effect of the electromagnetic potential due to a flux cf) together with the 
potential barrier supported on the infinite thin shielded solenoid. 

Our initial task is to provide a class of well defined operators corresponding to 
the formal differential (plus distributional) expression 

-j2(d/dx e -A(d e )) 2 -X6(r) , 

where x\, 22, ^3 are the standard coordinates in iR 3 , A = i(cj)/2iir 2 )( y — x-idx\ +Xidx^) 
is a pure gauge potential and r = ((a;i) 2 + (a^) 2 ) 1 ^ 2 . For that aim, we first reduce the 
problem to two dimensions by making use of translational symmetry with respect 
to the coordinate X3. Let r, 9 be, respectively, the radial and angular coordinate in 
]R 2 , < r and < 9 < 2n. Set a = —(J)/2tc. We concentrate on the case when 
a ^ 2L and owing to the gauge symmetry A' = A + e~ me de me , n 6 1, we can 
restrict ourselves to the case 

ae]0,l[. (1) 

The method we adopt in this paper is based on self adjoint extensions of sym- 
metric operators. From this perspective, we try to combine two well known cases 
which were already extensively discussed in the literature. The first one, with a = 0, 
corresponds to the so called point interaction (cf. [||]) and was studied in detail in 
|2|]. An operator in a one-parameter family is defined on a domain which is charac- 
terized by a linear relation between certain coefficients which are built up from the 
asymptotic behaviour, as r — > 0, of (singular) wave functions. The starting point 
was a symmetric operator with a domain formed by wave functions with supports 
separated from the origin. The deficiency indices turned out to be (1,1). 

The second case, with A = 0, corresponds to the pure Aharonov-Bohm (A-B) 
potential [l[ and was investigated many times on different levels (see 0, 0, ||], [§], 
|T0||). The generalized eigenf unctions are required to belong to H 2 ^(1R 2 \ {0}) and 



at the origin the regular condition 

lim/(r,0) = O (2) 

r— >0+ 

is imposed. As the Hamiltonian enjoys rotational symmetry the generalized eigen- 
functions are known and the scattering problem is solved explicitly. 

However, it has been known already for some time that when decomposing the 
Hilbert space into a direct sum, 

L 2 (M 2 , d 2 x) ~ L 2 (M + ,rdr) <g> L 2 ([0, 2tt], d9) = 0^ = _ oo L 2 (iR + , rdr) <g> Xm , (3) 

where 

Xm (9) = (2n)- 1 / 2 e md , (4) 
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then the A-B operator decomposes correspondingly and in the channels m = —1 
and m = the boundary condition (0) is not the most general one. But since the 
other boundary conditions admit wave functions which are singular at the origin 
they were usually ruled out. 

To our opinion, it makes good sense to consider the most general case and hence 
to allow even a sort of interaction between the two channels. Thus we apply to the 
pure A-B Hamiltonian exactly the same procedure which was used in the case of 
point interactions. The deficiency indices one obtains this way are (2,2) and this 
indicates clearly that the result is not simply a superposition of the two special 
cases. 

1 Five-parameter family of Hamilton operators 
and their resolvents 

In order to get operators which can be consistently interpreted as describing the 
physical situation we are interested in, we start with the case of pure A-B effect and 
introduce the point interaction at G IR 2 in the usual way. Namely, first we consider 
the restriction of the self-adjoint pure A-B operator H to the space of functions with 
supports outside of {0}, obtaining thus a closable symmetric operator. Then we shall 
find all possible self-adjoint extensions of its closure H. 

The adjoint H* is defined as the differential operator —(V — A(V)) 2 on the the 
domain 

^ G £>(#*) e L 2 (M 2 )nH^ c (M 2 \{0}) and (V — A(V)) 2 tp G L 2 (IR 2 ) . (5) 

On general grounds, H has equal deficiency indices. To find the corresponding 
deficiency spaces we employ the decomposition (||). Since the orthogonal projection 
onto L 2 (M + ) ®Xm commutes with H on T>(H), we can solve the eigenvalue problem 

H*f = k 2 f, (6) 

with k = e ,7r//4 and k = e* 37r//4 , in each sector m of the angular momentum. Setting 
f(r,9) = g(r)x m (0), (D becomes 

- (d 2 /dr 2 + 1/r d/dr + (m + af/r 2 )g = k 2 g, (7) 

which, by the standard substitution r — > kr, leads to the Bessel equation. 

Next, selecting in the two-dimensional space of solutions the one which vanishes 
at the infinity, we arrive at the Henkel functions 

9{r) = H^ +fA {kr) . (8) 

To ensure the integrability we still have to control the asymptotics as r — ► 0+. The 
case a = is known; the L 2 solutions exist only in the sector m = and thus the 
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the deficiency indices are (1,1). Assuming now that < <fi < 2tt and recalling the 
asymptotics 

HU(z) = + . t, ^ + 0(z-^) , (9) 

sin 7ns 1 (1 — jy) sin iru 1 {\ + v) 

as z — > 0, the integrability at means that 2|m + a| — 1 < 1, which selects precisely 
two angular momentum sectors: m — — 1 and m — 0. Thus the deficiency indices of 
H are (2, 2) and the deficiency space Mi is spanned by f} and given by 

/, l (r, (9) = NiH&We* , /, 2 (r, (9) = N 2 H^(kr) , (10) 

where = e l7r//4 (fc = y/i with Im > 0) and the normalization constants Ni,N 2 
will be determined later on. 

This means that all self-adjoint extensions are in one-to-one correspondence with 
the elements of the unitary group U(2) and are determined by boundary conditions 
at the origin. We treat them in detail in the next section. Thus we get, apart of a 
characterizing the magnetic flux, four additional parameters. 

It is now straightforward to determine the domain T>(H). As H = H** it holds 
true that 

if) G "D(H) -<=>• ip G D(H*) and (ip, H*(p) = (H*ip,(p), Vtp G K + M- t . 

Consequently we find that g(r) <8> Xm{9) G T>(H), with m = —1, 0, if and only if 
g G L 2 (R + , rdr) n fl^(]0, +oo[), (<9 r 2 + r" 1 ^, + r~ 2 (m + a) 2 )g G L 2 (iR + , rdr), and 

Urn irW{g,h±) = 0, (11) 

r— >0+ 

where h± = H$ +Ct ,(y/5r) (with Im y/S > 0) and the symbol W(g, h) stands for 
the Wronskian, 

W(g, h) := gd r h — hd r g . 

At this point, let us make a short digression and recall a useful formula comparing 
resolvents of two self-adjoint extensions. It is stated in the framework of Krein's 
approach to the theory of self-adjoint extensions and was presented originally in |J. 
Let us consider a general situation when a closed symmetric operator X is given and 
Aq is a self-adjoint extension of X, X C Aq = A* Q C X* . Assume that the deficiency 
indices of X are {d, d), d < oo. Set 

M z ■= Ker(A* - z), R° z := (A Q - ze€\M. (12) 

The following facts are well known and easy to check. First, 

i> eM w i> + (z - w)R° z i) = (A - w)R° z i) G M z ■ (13) 

Second, if U® : M z — > Mz, z G(C\iR, is the unitary mapping defining the self-adjoint 
extension A then 

V^GA4, U° z ^ = -(A -z)(A -z)- 1 ^. (14) 
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Fix w gC \ M and a basis {f w }t in Af w . Set 

f e z = fi + (z-w)FP x fi = (Ao-w)SP g fi, zeC\M. (15) 

Since (A — w)R° z is injective and dim Af w = dimA/" 2 = d, {f z , . . . , f z } is a basis in 
Af z . Suppressing the superscript £ one can verify readily that 

Vz, z gC \ M, f z = f z , + (z-z')R° z f z ,, (16) 
VzeC\lR, U°J z = -f- z . (17) 

Thus, in order to reproduce the vector-valued function f z , one can take any z' GC\iR 
instead of w. Furthermore, the matrix of U z in the bases {f z }e, {fg}e equals — 1. 

To proceed further let us introduce a matrix P(z,w) = (P^ k (z,w)) of scalar 
products relating two spaces Af z and Af w , 

P jk (z,w):=(fi,f*). (18) 

One finds that 

P(z,z)=P(z,z), (19) 

and if U is any matrix expressed in the bases {f z }, {/!}, and representing a unitary 
mapping Af z — > Af z then 

P(z,z) =U*P(z,z)U . (20) 
Furthermore, if P z is the orthogonal projector onto Af z then 

p 2 =E(^^) _i yVi(/^-)- ( 2i ) 

Our primary interest is to compare A with another self-adjoint extension A of X 
corresponding to a family of unitary mappings U z : Af z — > A/" 2 . Krein's formula tells 
us that 

R z -R° z = (z- z)- x V* z {U- z - U° Z )V- Z , (22) 

with the symbol V* standing for the embedding of M z into the global Hilbert space. 
This means that there exists a family of d x d matrices, M z = (M| fc ), such that 

R z -R z = Y, MJ z k fi </£•>■ (23) 

We claim that 

.U. - .U„. (: - w)M z P(z, w) M w , (24) 

(iu-«;)M w = (C/- 1 + l)P(«;,«;)- 1 , (25) 

where U is the matrix of U w : A4, — > Afw in the bases {/^} and {f w }. The proof is 
quite straightforward and relies on Krein's formula, the first resolvent identity and 
the explicit expression for V z . Provided M w and M z are invertible one can also write 

M- 1 -M- 1 = (z-w)P(z,w). (26) 
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What we shall need in the sequel is the particular choice of w = —i. In this case, 

M" 1 = 2iP(i, i)(U + I)" 1 -(z + i)P(z, -i) , (27) 

where this time U is the matrix of U % : J\T% — > M- % in the above specified bases. 

Next we apply this general procedure to our problem, with A = H - the pure 
A-B operator. Thus, for z G C \ M + , we choose in M z a particular basis which 
depends holomorphically on z by 

fl = {l + {z-%)R M )fl £ = 1,2, (28) 

where i?^ is the resolvent of the pure A-B operator defined by its integral kernel 
(Green function) 

G z (r,e;r',e') = — £ ^ } / {k 2 - z)- l J\ m ^{kr)J ]m + a \{kr')kdk . (29) 
Recalling that 

HW(z) = -e-™'/ 2 K v {-iz), (30) 
m 

where 

K v {z)=\ e~ zcosht coshut dt, (31) 
j o 

and using the formulae 8.13(3) and 8.5(12) of 

^ if M (oy)^(«y)ydy=^-J (z' + a 2 )- 1 . ( 32 ) 



/ -= -J t ,(yx)xdx = a>*K lx (ay), /i e } - 1, 3/2 [, Re a > 0, (33) 



we have 



£(r,0) = iV ie — ( 1 - Q )/ 4 ( v ^) 1 - a i7 1 a) Q (v / ^r)e- ie , 

/ z 2 ( r ,#) = N 2 e- ma /*(y/z) a HgXVZr) , (34) 

where z 6(C\ iR + , Im y'i > 0. Particularly, 

/i (r, 0) = N^-^^H^kry- 19 , /^(r, 0) = N 2 e ma/2 {kr) , (35) 



where /c = = e l3n ^. 

A word of warning should be said here. We use the branch (e tip ) u = e lipu , for 



< (p < 2tt, which differs from the usual choice made in surveys like ||, |fLl 
Next we compute the matrix P(z, z'), 

P e ' m (z, z') := {f z , f?) = P e >\z, z')5 e>m , £,m = l,2. (36) 



Since K v {z) = K v (z) and using the formula 6.521 of 0, 



K„(ar)K 1/ (br)rdr = \ > \ > 37 

2 sm(7rz/)(cr — tr) 
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we find that 

(38, 

sm(7n/J 2! — 2 

where 

_ J 1 - a if £ = 1 
U ~\a if £ = 2 ■ 

Making use of the identity 

J \Hjp(kr)\ 2 rdr = (tt cos(7rz//2)J , fc = \/5, Im > 0, (39) 
we choose now 

iVi = 2- 1/2 sin 1/2 (7ra/2), N 2 = 2~ 1/2 cos l/2 (na/2) . (40) 

In this case the basis {//} (|10|), as well as the basis is orthonormal, i.e., in 

the matrix notation, 

P(z,z) = P(-t,-«) = J , (41) 
where 7 is the 2x2 unit matrix. Moreover, introducing the matrix 

D = diag(l — a, a) , (42) 

we also have 

-i) = (z + z)" 1 sm- 1 (nD/2)e- mD (e 3mD/2 - z D ) . (43) 

We conclude that all self-adjoint extensions H u of H are bijectively labelled by 2 x 2 
unitary matrices U by means of 

V(H u ):={t, ^ = "0 + £q^}, ( 44 ) 

where ?/> G D(H) and 

m 

and 

# ^ = ^ + E - » E /-^ m,i ) • (45) 

t m 

According to the above discussion, U = —1 corresponds to the pure A-B operator 
H. Moreover, diagonal U describe the extensions preserving the angular momentum 
(which otherwise is non-conserved). 
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2 Boundary conditions 



The family of operators H u , defined so far abstractly, can be equivalently charac- 
terized as differential operators with some well specified boundary conditions. For 
this purpose, we introduce four linear functionals n — 1, 2, a — 1, 2, correspond- 
ing to the two critical angular sectors and to the first two leading terms giving the 
asymptotic behaviour of the radial part of ip as r — ► 0. We define 



*i(V0 


:= lim r 




:= lim r 

r-»Q 


*?(V0 


:= lim r 


*2(V0 


:= lim r 

r->0 



1-a 



2tt 



%(r,#)e* e ^/2tt, 

2 " ip(r, 9) e ld d9/2n - r" 1+a $}(^) 
^(r,6)de/2n, 
2W ip(r,9) dB /2tt - r~ a $?(V>) 



(46) 



This definition is, of course, dictated by the asymptotic behaviour of Hankels func- 
tions (cf. (^)). So for t/> G V(H u ), the part of ■?/> which is singular or becomes 
singular after differentiation by <9 r is given by 

Let us first check the symmetry condition (<£>i, H U {p 2 ) = (H U (pi, <p 2 ), for <£>i, G 
T>(H U ). The integration by parts gives (W = Wronskian) 



lim 

r^0+ Jo 



2 7! 



rW((p U (p 2 )dO = 0. 
Only the singular parts of (fx, (p 2 can contribute and thus one arrives at 

$i(pi)*D$ 2 M = V^ 1; ^ 2 e V(H U ) , 

where we have introduced 

iW := 



1,2, 



(47) 



(48) 



and D was defined in 

Next we apply the functionals to the functions Namely, introduce four 
matrices $ a 6, where the label a = 1, 2 refers to the first (respectively the second) 
leading coefficient, and b — ± refers to ± in A/± l . They are defined by 



(49) 



so the rows of these matrices are numbered by the angular momentum n (for the 
sake of convenience we shifted the index by 2, n = m + 2G {1,2}) and the columns 
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by £, corresponding to the basis in A/± 4 - In view of the asymptotic expansion of the 
functions (cf. ([[]]) and (pop), they read explicitly 

= — = — - cos 1/2 (tt j D/2) 2 d T(l - D)- 1 exp(-t7cD/A) , 
sin(jra) 

oO- 1 / 2 

$2+ = — 7 — r cos 1/2 (7r J D/2)2~ D r(l + J D)- 1 exp(-i37rZ>/4) , (50) 
sm(7ra) 

$ 2 _ = — 7 — r cos 1/2 (tt j D/2) 2~ d T(l + D)- 1 exp(inD/4) . 
sm(7ra) 

Here and everywhere in what follows we use the obvious notation: if a function F 
is well defined on the set {1 — a, a} then 

F(D) := diag(F(l - a), F(a)) . 

From the definition (0) of V(H U ) follows that ip E L 2 (R 2 ) belongs to the 
domain of H u if and only if 

Inversely, assume that we are given a couple of 2 x 2 matrices X\ , X2 such that 
rank(X',X2) = 2, and consider the boundary condition 

^ERan^'V (52) 



*2«0/ VA' 2 
The symmetry property ( 47|) leads to the requirement 

X*DX 2 = X*DX! . (53) 

In fact, relying on the explicit form of the matrices $ a f,, one can show quite straight- 
forwardly that for any couple Xi,X 2 with the above properties there exist exactly 
one 2x2 unitary matrix U and Y G GL(2,€) such that 

X X Y = + x 2 r = $ 2 ,+ + $ 2j -^ • (54) 



On the contrary, if U is unitary then X a = $ a + + $ a _U, a = 1,2, verify ( p3| ) 
and rank(X',X2) = 2. This way we have rederived a well known result that all 
self-adjoint extensions of H are in one-to-one correspondence with points of a real 
4-dimensional submanifold of the Grassmann manifold G^fC 4 ) determined by the 
equation (|53|). 

One can rewrite the boundary condition (|52"D in a more convenient form when 
making use of the biholomorphic diffeomorphism G^C 4 ) - ► G^^C 4 )"- Here G^C 4 )" 
stands for the Grassmann manifold in the space dual to C 4 . The points of G^lp 4 )" 
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are represented by couples of matrices A\,A 2 G Mat (2,C), with rank(A 1 , A 2 ) = 2, 
modulo the left action of GL(2,(C). The diffeomorphism is given by 

Ran ) e G< 2 (€ 4 ) i- Ran ) G G 2 (C 4 ) , 

where A X X X + A 2 X 2 = . (55) 

The real submanifold of G 2 (([7 4 ) determined by (|53| ) is mapped bijectively onto the 
real 4-dimensional submanifold of G 2 (([7 4 )^ determined by 

A 1 D- 1 A* 2 = A 2 D- l A\. (56) 

We conclude that each self-adjoint extension of H is determined by a boundary 
condition of the type 

A 1 $ 1 (^) + A 2 $ 2 (^) = 0, (57) 

where A U A 2 G Mat (2,C) verify (||) and rank(Ai,A 2 ) = 2. Two couples (A 1 ,A 2 ) 
and (A^,^) define the same self-adjoint operator if and only if there exists Y G 
GL(2,€) such that A[ = YA U A' 2 = YA 2 . 

Let us now restrict ourselves to an open dense subset of the manifold of boundary 
conditions (|57|) which we obtain by fixing A\ = I and setting A 2 = —A. The 
condition ( |57| ) means, of course, that 

The restriction ( |56"D then reads 

DA = A*L> . (59) 

All matrices A obeying (|59| ) can be parameterized by the aid of four real (or two 
real and one complex) parameters, namely 

A = ( U , aW ] , with u,v e R and w e@. (60) 
V(l — a)w f / 



From (|55| ) follows that one can choose Xi = A, X 2 = /. In virtue of (|5^), this 
leads to a relation between A and U, 

A = ($i,+ + $ li _E/)($ 2 ,+ + $2,-f/) _1 , (61) 

provided the relevant matrix is invertible. We use the following parameterization of 
U, 

ge ta _( 1 _ g 2^ ) l/2 e - l 6- 

where a, b, q, to G M, < g < 1. The matrix ($ 2 ,+ + $ 2i _?7) is invertible exactly 
when d ^ where 

<i := sintu + gsin(a — na) . (63) 



" = W 1^ I- (62) 



10 



In this case, the parameters of the matrix A can be expressed explicitly in terms of 
a,b,q,u, 



-1 o2-2a r ( 2 -«) / / . , K\ , _ / 7T 



u = d 2 a — lcos[uj + —aJ+qcos\a a 



cT 1 (2(l-g 2 )sin(vra)) 1/2 e-^^), (64) 



d- 1 (2(l-g 2 )sin(7ra)) 1/2 e< & -5^), 



i-1 2ry r(l + a) / / 71 \ / 71 

v = — d 2 — sin w a +osin a a 

r(l-a) V V 2 / H \ 2 

Obviously, w = w = v = corresponds to the pure A-B case. Moreover, diagonal 
A describe the operators preserving the angular momentum (w is responsible for its 
non-conservation) . 

3 Spectrum and eigenspaces 

In order to find the spectrum one can use Krein's formula for the resolvent B% = 

R U g =R z + J2f^Mf(fi,.) . (65) 

kt 



Using (0), (H) and (g7|) we get 

M~ l = sin- 1 (7r£>/2) (VV i7rD ([7 + 1) - e lnD/2 U - e -' mD/2 ) (U + l)" 1 . (66) 

Since R v z is a rank two perturbation of R z , H u and H have the same absolutely 
continuous spectrum, namely [0, oo[. 

The discrete spectrum is determined by the equation det M~ l = 0, i.e., 

det (p 2D (U + 1) - e mD/2 U - e - lnD/2 ) = , (67) 

where we have introduced p > by p 2 = —z. According to the discussion below 
there are no non-neg ative eigenvalues. If d ^ (cf. ©) then (e mD ' 2 U - e~ mD ' 2 ) 
is invertible and owing to (pTj), 

(U+ \){e mD / 2 U -e-^ 012 )' 1 = - cos- 1 / 2 (7rD/2)2- D r(l- J D)e 47rD/4 A 

x cos 1/2 (7rD/2)2- D r(l+D)e^ 7rD/4 . (68) 

In this case, (|67|) is equivalent to 

det fEMW !D A + 1 ).„. ( 69 ) 



T(l+D) V2 
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Rewriting fl69]) in terms of the parameters u,v,w (cf. (|60|)) we get 

~ 2a r(i-a) ^ 




rri+a 



w = |wr. (70) 



Consider the LHS of ( |70|) as a function defined on ]0, +00 [. Since the RHS of 
(|70|) is always non-negative some elementary analysis gives immediately the number 
of solutions. The number of roots of F(p) equals or 1 or 2. There is no root iff 
u > 0,v > 0, and there are two roots iff u < 0, v < (it may happen that the 
two roots coincide giving a multiple root). Denote by p\ the smallest root of F(p), 
if any. In the case of two roots, let P2 be the greater one. Then clearly F(p) is 
decreasing on ]0,pi] from +00 to zero and is increasing on [p 2 , +oo[ from zero to the 
asymptotic value uv/a(l — a). 
We conclude that there are 

two eigenvalues if u < 0, v < 0, and det A = uv — a (1 - a)|w| 2 > 0, 
no eigenvalue if u > 0, v > 0, and det A = uv — a(l — a)|w| 2 > , (71) 
one eigenvalue otherwise . 

We stress once more that all eigenvalues, if any, are negative. Again, it may happen 
(when w — 0) that two eigenvalues coincide producing consequently a multiple 
eigenvalue. 

For generic a not much can be said about what are the solutions, except the 
case w = when p = 2 (~ r(2^) ) if w < and/or p = 2 (— Yct^f ) if t> < 0, 
and the case u = v = when p = 2/\w\. An interesting particular case is a = 1/2 
when we can give a complete answer about the values of the two solutions (the first 
case in (|7lD ), 

_ u + v±(\w\ 2 + (u-v) 2 y* 

p± ~ 2(M74-H ' 1 ' 

and of the solution (the third case in ([71])), 
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= u+v+(\w\ 2 + (u- , 
P 2{\w\ 2 l±-uv) ' K ' 

Similar, but more complicated analysis can be also performed for a = 1/3, 1/4 and 
partially for other fractional values of a. 

As far as the eigenvectors are concerned, they have to be obtained, of course, 
as solutions of the differential equation (D) including the corresponding boundary 
conditions Q58|). First of all, it is clear that in the sectors of the angular momentum 



m 7^ —1,0, there is a complete system of generalized (and normalized) eigenfunctions 
coinciding with those of the pure A-B effect, 

(2n)- 1 / 2 J lmHlll (kr)e me , m = . . . - 3, -2, 1, 2 . . . , k > . (74) 

Next we pass to the sectors m = — 1 and m = 0. 
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As far as the (true) eigenfunctions are concerned, the L 2 -integrability condition 
at infinity restricts the eigenvalue k 2 to k 2 < 0, and picks up a unique solution, up 
to a multiplicative constant, in each sector m = — 1 and m = (this means that 
both exponential growth and oscillatory behaviour at infinity are excluded). Hence 
setting as before k = ip, with p > 0, the eigenfunction must have the form 



iH { ^} a {ipr)e- lt) + r)H^{ipr) 



(75) 

where G W. The boundary conditions (|58D lead to the following relation between 
i and 77, 



(l+{^) D T{l-D)KT(l+D)^e-™ D (f 



. 



(76) 



Setting the determinant of this system of linear equations to zero we recover the 
equation on eigenvalues (|70|) . Any non-trivial solution (£,r/) of ([[6|) determines an 
eigenfunction in accordance with ([75]). 

As far as the generalized eigenfunctions are concerned, the eigenvalue equation 
admits a four-parameter solution 



{£J-i+ a (kr) + v-Ji-aM) e' ie + e-J-a(kr) + V f J a (kr) 



(77) 



where k > and £, r/, r/' G (T. 

In view of the asymptotics of J u (z) ~ r(l + z/) _1 (z/2)^ (1 + 0(^ 2 )) and again by 
taking into account the boundary conditions we find the relation 



A possible choice is 



{k/2) D T{l-D)AT(l + D)- 1 (k/2) 



d ( V 
V' 



(78) 



and 



and in this way we obtain two independent solutions 



where 



b 2 (k) 



faJ-i+aikr) + Ji- a (kr)) e~ %e + Vl J_ a (kr) 
&J-i+ a (kr)e- ld + T] 2 J^ a (kr) + J a (kr) , 



(79) 



(k/2) D T(l-D)AT(l+D)- 1 (k/2) 

2-2a 



D 



( V{a 

u 



r(3-a) V2 
k 



w ■ 



T(l-a) 



k 



w 



Set 



N(k) 



6 6 
»7i m 



+ exp(ziiD) . 



2a 



W) 



51) 
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Now we seek a pair of eigenf unctions which are complete and orthonormal in the 
generalized sense. In order to compute the scalar products of b\{k) and 62(^)5 we 
need to know the integrals involving the products J^ay) J^{xy) and J_ tl {ay) J^xy). 
Recalling the relation between the functions H^' and fl3"0|) and using the limit 
value of (|32T) we have 



DC 



H V (ay). Uxy)ydy = — (-)' ( l — - -J- ) . (82) 

^ ma \aj \x — a — tu x+aj 

Next, with the help of the distributional identity (x — iO)' 1 = V(l/x) + i7r8(x), 
where V denotes the principal part, we get two identities by taking the real and 
imaginary parts of (p2|) , 

J^{ay) J il (xy)ydy = -S(x-a), (83) 

J-^ay) J fl (xy)ydy = C ° S ^^ 5(x-a) 

sin(7ru) 1 \ . . 

+ — — - \V{ • 84 

ixa \aj \ x — a x+aj 

Applying ( |33|) -(|84|) to the solutions (|79| ) arranged in a row B(k) = [bi(k), b 2 {k)), 
we obtain the following 2x2 matrix of scalar products, 

(B(k'), B(k)) = N(k')*N(k)rS(k - k') . (85) 

k 

We observe that 



uv 



, A . , . ... \ k 2 uT(a)e ma k 2 - 2a v T(l-a) e~ ma tf a , 

det Nik) = — r — ww I — — h „\ ' , r \ — '- r 1 86 

v ; \a(l-a) J 4 2 2 ~ 2a F(2-a) 2 2a T(l+a) V ; 

equals minus the LHS of (|69|) , with p being replaced by —ik, and thus, in view 
of our analysis of eq. (0), det N(k) is nonvanishing for all k > 0. Therefore 
(g 1 (k),g 2 {k)) := B(k)N(k)~ l , being given by 

gi (k;r,e) = det" 1 iV(A ; )((^2-6^ + eie i7ra )J a -i(A;r)e-^ 

+ (V2 + e™ a )Ji-<*(kr)e-* 6 - ^J^kr) + Vl e ma J- a (kr)) , 

g 2 (k; r, 9) = ter x N{k) {-^'^J^kr^ 6 - kJ^{kr)e* (87) 
+ (6 - e~ ma ) J a (kr) + ~ km ~ V2e- l7Ta )J~a(kr)) , 



where £1, £2, Vi an d V2 are defined by (p0|), form a complete orthonormal basis of 
generalized eigenvectors in the subspace corresponding to the absolutely continuous 
spectrum in the two considered sectors, 

(g j (k'),g t (k)) = ^S{k-k')S j4 . (88) 
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4 Scattering 



The existence of a complete and orthonormal basis of generalized eigenvectors is 
sufficient to show that the wave (M0ller) operators W± = lim^-i-oo e ltH e - ltH ° exist 
and are complete. In fact, they can be exhibited explicitly as well as the scattering 
operator S = (W + )*W_. 

From what we have said so far it is evident that W± and S preserve the sectors 
m 7^ —1, 0, and there, they are exactly the same as in the pure A-B case; in particular 
S m = e 2l<5m , where 5 m = (|m| — \m + a\)w/2. Thus we restrict ourselves to the 
subspace L 2 (M + ,rdr) <g> (Cx-i © <Cx-o) > of remaining sectors m = —1,0, which 
is, of course, also preserved by all the relevant operators (if there is no danger of 
confusion we denote the restriction of an operator by the same symbol). 

Using the basis g±(k) , g 2 (k) we define a unitary operator T from L 2 (M + , kdk)&C 2 
to L 2 (R + ,rdr) ® (fx-i ©<Cxo), by 

_ ^ poo _ 

^]=V:=E/ gAWjiQkdk, (89) 

with the inverse J 7 " 1 [■?/>] = ip given by 

^(k) := (g 3 (k)^) . (90) 

The operator T satisfies 

e- ltu T\i)\ = T\e~ lt ^ 4>) (91) 

We append the superscript t0 ' to the relevant objects like g^{k) or JF° corresponding 
to the free case (a = 0, A = 0); particularly 

9x(k; r, 9) = Ji(Ar) e~ % \ g° 2 (k; r, 9) = J (kr) . (92) 

Now we seek a pair of 2 x 2 matrices Q± = Q±(k), generally depending on k and 
acting in an obvious way as a multiplication operator on L 2 (M + , kdk) ®C 2 , so that 

W± F° = T Q± . (93) 

Then it follows that S J 1 " = JF° and this means that X := is nothing 

but the scattering matrix in the momentum representation (restricted to the sectors 
m = -1,0). 

We have to verify that 

lim \\e- ltH ° J*M] - e- ltH F\n±ijA\\ = . (94) 

t — >±oo 

Due to ([H]), the condition fl9~4]) means that 

poo „ 

lim / ^ \(9-(k) -Y,9l(k)n ± (k) tj )1>3(k)\ 2 k dk = °- ( 95 ) 
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It is sufficient to prove fl9"5|) only for functions ipj{k) from the dense subspace 
C^°(M + ) C L 2 (M + , k dk) . By the stationary phase method, the convergence of 
such an integral, as t — > ±00, will be established provided the coefficient standing 
in front of the term e ±lkr vanishes. In view of the known large x expansion 

J„{x) = (27rx) 1/2 (e i(:E -^ /2 -" /4) + e -*(*-*W2-ir/4)) + 0(x 3/4 ) , (96) 

we obtain, by looking separately at the coefficients in front of e~ %e and 1, that 

^(A;)- 1 = diag(-e- mQ/2 , e ma/2 )N(k)N(ky l 
= det' 1 N{k) diag(-e~ ma/2 , e ma/2 ) 
/fir?9 — £9771 



£1^2 - - e 4 ™(r/ 2 - 6) - e 2ma (e mQ - e" ma )^ 



2 

(e mQ - e- 4 ™)//! ^772 - ^ - e- ma (r] 2 - ft) - e- 2 *™ ,/ ' 

(97) 

where 

^(*)=(^ ^J+expHvrD), (98) 

and also 

fi_(ife) = diag(-e- 4 ™ /2 , e J7ra/2 ) . (99) 

Note that 

N{k)*N{k) = N{k)*N{k) (100) 

and so Q + (k) is unitary. Consequently, we can express (the entries of) £ in terms 
of the parameters u,v,w (see (|S0D) as 

S n = det _1 iV(fc) x 

( P -™ ( UV - \ W A k l + uV ^)k 2 - 2a _ vT(l-a)k 2 « _ \ 
\ \a(l-a) 1 1 J 4 2 2 " 2Q r(2-a) 2 2 «r(l+a) J 

E12 = — det _1 iV(/c) 2 sin(7ra) wk 

E 21 = -det _1 iV(£;) 2 sin(vra) wfc (101) 
S 22 = det~ 1 A^(A;) x 



W uv 2 U 2 ^(a)^ ^(l-q)fc 2 " _ ;jr „ 

' 1 a(l-a) 1 1 I 4 2 2 - 2a r(2-a) 2MYl + ai) ' 



where det _1 iV(fc) is given by (|S6D. Obviously, £ is unitary. 

Let us specialize these formulae to three particular cases, 
(i) If w = (conserved angular momentum) then 



£ = diag 



'uT(a)(k/2) 2 - 2a - e ma T(2-a) vT(l-a){k/2) 2a + e~ ma T(l + a) 



e ma uT(a){k/2) 2 - 2a -T(2-aY e- ma vT(l-a)(k/2) 2a + T(l+a) ) ' 

(102) 

(ii) If u = v = (maximal nonconservation of angular momentum) then 

y _ n -i (e- ma \w\ 2 {k 2 /4) + e ma ism{-Ka) wk \ , , 

Q X { isin(na)wk e wa \w\ 2 {k 2 /4) + e~ ma J ' [ ' 
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where 

Q = \w\ 2 ^- + l. 

1 1 4 

(iii) If a = 1/2 then 

£, _ „_! / — 2 — i(uv — \w\ 2 /4)k 2 + (tt — t>)A; — mD/c \ 

V — m>£; z + i(uv — |w| 2 /4)/c 2 + (u — v)k J ' 

(104) 

where 

Q = -l + (uv- \w\ 2 /4)k 2 + i(u + v)k . 

We conclude this section by giving (the kernel of) the full scattering operator in 
the angular representation, 

1 oo 

S(k-9,9 ) = — £ S(k;m,n)e< me - n ^ , (105) 



m,n=—oo 



where 

e * W (|m|-|m+a|) ^ jfm^^-^O 

S m +2,n+2 ifra,ra = -l,0 



S^fc; m, n) = 



(the shift by 2 is due to our labelling of rows and columns of S). The double infinite 
sum can be performed in the sectors m,n ^ —1,0 (a known result borrowed from 
the pure A-B effect) and yields 

X ^ e -i(0-0 o )/2 

S(k;9,9 ) = cos(7ra)5(9 — 9 ) H sin(7ro; N 



2tt v sin(fl/2 - o /2) 
+ J- ^ (E m+2 , n+2 - e-( 2m+1 >™5 mn ) e »(""-"*0 . (106) 



m,n=— 1 



We recall that the differential cross-section in the plane is given by da(9)/d9 
(2n/k)\S(k;9,9 )\ 2 . 



5 Conclusions 

We have introduced and studied a five-parameter family of operators which describe 
a quantum mechanical particle interacting with a magnetic flux a caused by an in- 
finite thin material solenoid. We conjecture that in some well defined limit (when 
the thickness of the solenoid — > and its length — > oo) which should be largely 
independent on the details of approximating potentials (and which however goes 
beyond the scope of this paper) such a situation is described by a (singular) poten- 
tial barrier and by a electromagnetic potential concentrated along the z axis (the 
magnetic field vanishes in the remaining region). One of the five parameters is just 
the value of the flux a and the other four correspond to the strength of a singular 
potential barrier (sort of a combination of Dirac 5 and 5') and can be interpreted as 
penetrability coefficients of the shielded solenoid. 
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A general operator of this family corresponds to an intricate mixture between 
the Aharonov-Bohm effect and the point interactions which is manifested more con- 
cretely via the mixing between the angular and the radial boundary conditions. It 
is interesting that the result we have obtained is richer than a simple superposition 
of the point interaction and the pure A-B effect. For instance, for a range of pa- 
rameters there are two bound states possible while for the usual interactions with a 
support concentrated along the z axis and symmetric under the ^-translations (or, 
equivalently, for a point interaction in two dimensions) there is always (except of the 
free case) exactly one bound state, and for the pure A-B effect there are no bound 
states at all. 

In the present paper we have derived an explicit formula for the scattering matrix 
S(k; 9, 9q) depending on the five parameters. Naturally, it would be of interest to 
examine the differential cross-section numerically, particularly its dependance on 
the parameters, and to deduce some physical consequences. This is what we plan 
to do separately. 

It is a matter of experimental measurements (interference, scattering or in con- 
densed matter) to establish which of the extensions in our family correspond to 
realizable models. On general grounds, one can distinguish some class of exten- 
sions by eg. symmetry requirements (conservation of the angular momentum), or 
postulating that there are no bound states (solenoid as a repulsive barrier). 

After completing most of our work we became aware of which has an overlap 
with some of our results. We have been also informed that a related preprint by R. 
Adami and A. Teta is in preparation. 

Acknowledgements. P.S. wishes to gratefully acknowledge the partial support 
from Grant No. 201/94/0708 of Czech GA. 
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